Radio Galaxy Clustering at z ~ 0.3 
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ABSTRACT 

Radio galaxies are uniquely useful as probes of large-scale structure as their uniform identi- 
fication with giant elliptical galaxies out to high redshift means that the evolution of their bias 
factor can be predicted. As the initial stage in a project to study large-scale structure with 
radio galaxies we have performed a small redshift survey, selecting 29 radio galaxies in the range 
0.19 < z < 0.45 from a contiguous 40 deg 2 area of sky. We detect significant clustering within 
this sample. The amplitude of the two-point correlation function we measure is consistent with 
no evolution from the local (z < 0.1) value. This is as expected in a model in which radio galaxy 
hosts form at high redshift and thereafter obey a continuity equation, although the signahnoise 
of the detection is too low to rule out other models. Larger surveys out to z ~ 1 should reveal 
the structures of superclusters at intermediate redshifts and strongly constrain models for the 
evolution of large-scale structure. 

Subject headings: large-scale structure of the Universe — galaxies: active — surveys 



1. Introduction 

Powerful radio sources are almost exclusively 
associated with giant elliptical galaxies, and ap- 
pear to be in richer than average environments 
(e.g. Hill & Lilly 1991). This suggests that they 
should be more biased tracers of the mass distri- 
bution than normal galaxies. A study of the clus- 
tering of local (z < 0.1) radio galaxies by Pea- 
cock & Nicholson (1991) showed that this was in- 
deed the case, with radio galaxies having a cross 
correlation function of the usual form (assumed 
throughout this paper), £ g9 = (r/i"o) , with a 
correlation length of ro = ll/i -1 Mpc. 1 This 
can be compared to ro = 5.7ft. Mpc for normal 
galaxies (Loveday et al. 1992); ro = 4.5/i~ 1 for 
IRAS-selected galaxies (Fisher et al. 1994), and 
ro = 14.3/i _1 Mpc for rich clusters of galaxies (Dal- 
ton et al. 1994). Radio galaxies thus cluster with 
a strength intermediate between normal galaxies 
and clusters (Bahcall & Chokshi 1992). The as- 



1 We assume a cosmology with 
100/i _1 kms~ 1 Mpc~ 1 throughout 
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sociated cosmological bias factor is lower than for 
clusters, but about twice as high as that for IRAS- 
selected galaxies (Peacock & Dodds 1994). 

Studies of the angular correlation function of 
radio sources by Magliocchetti et al. (1999) and 
Cress & Kamionkowski (1998) show that the data 
from radio surveys such as FIRST are consistent 
with little evolution in the clustering amplitude 
in the range < z < 1. However, these conclu- 
sions are based on extrapolating the luminosity 
functions of Dunlop & Peacock (1990; hereafter 
DP), to the faint flux densities near the limit of 
the FIRST survey, S1.4 ~ 3 mJy. At these levels 
the DP luminosity functions are constrained by 
source count data only, as direct redshift surveys 
were only available at Si. 4 > 200 mJy. In particu- 
lar, at faint flux levels the radio source population 
is a mix of nearby star-forming galaxies and AGN- 
powered radio sources with a range in redshifts 
from to >4 (Condon et al. 1998), which have 
quite different clustering properties. It is therefore 
important to test the results of the angular corre- 
lation function studies with direct measurements 
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of clustering from radio galaxy redshift surveys. 

Studies of the clustering of radio-quiet AGN 
seem to show a generally similar correlation 
length, but there is a wide range in estimates 
of r from different samples. This can probably 
be explained if the correlation function depends 
both on redshift and AGN luminosity [e.g. Sabbey 
ct al. (2000), La Franca, Andreani & Christiani 
(1998)]. 

Magliocchetti et al. (1999) discuss theoretical 
predictions for the evolution of the two-point cor- 
relation function of radio sources. Perhaps the 
most appropriate case to take is that where ra- 
dio galaxies form at high redshift (z 3> 1). We 
can trace the evolution of the host population 
out to z ~ 3, and find that the hosts vary lit- 
tle with redshift, apart from some passive evolu- 
tion. The host magnitudes are also only weakly 
dependent on radio luminosity (Lacy, Bunker & 
Ridgway 2000). Hence uncertainties in the evolu- 
tion in the bias factor are unlikely to be as im- 
portant an issue for radio galaxies as they are for 
normal galaxies or radio-quiet quasars. Fry (1996) 
shows that in this "galaxy conservation" scenario, 
the bias factor increases with redshift according 
to b(z) = 1 + (bo — 1)(1 + z), where bo is the bias 
factor at the present epoch. This is because fluc- 
tuations in the galaxy density field are fixed at 
the epoch of formation, but the fluctuations in the 
matter density field grow with time. The decrease 
in the bias factor with time is mostly compensated 
for by the clustering of matter under gravity, for 
which the growth factor D(z) = (1 + z)^ 1 for an 
CIm — 1, = cosmology. The two-point cor- 
relation function, cx D 2 (z)b 2 (z), should therefore 
show little evolution. 

To test this model, and to examine the nature 
of intermediate redshift superclusters, we therefore 
decided to begin a survey of large-scale structure 
at moderate redshifts (z <~ 0.2 — 0.65). This paper 
describes the initial result from this survey and 
also the prospects for future surveys. 

2. Survey strategy and observations 

We tried to optimize the survey to detect 
supercluster-scale objects at z ~ 0.4. We there- 
fore picked a point on the radio luminosity func- 
tion where the space density of objects close to 
the flux limit would be ~ 10~ 5 /i 3 Mpc~ 3 , thus ob- 



taining several objects in structures of linear sizes 
~ 100/i _1 Mpc. This corresponds to a radio lu- 
minosity of « 2 x 10 23 h~ 2 WHz _1 sr _1 , or a flux 
limit of « 20mJy at 1.4 GHz. This is comfortably 
above the completeness limits of the FIRST and 
NVSS radio surveys (Becker et al. 1995; Condon 
ct al. 1998). 

Initial selection was made from the NVSS cata- 
logue, each NVSS source was examined in FIRST 
to check for confusion, and to estimate the position 
of the identification. This technique combines the 
sensitivity to extended flux of the NVSS survey 
with the positional accuracy of FIRST. The sam- 
ple discussed in this paper consists of 322 objects 
within right ascension and declination ranges of 
01 h < R.A. < 01 h 48 m and -02° < Dec. < 01°20' 
(an area of ~ 40deg 2 ). 

Identifications were made on the UK Schmidt 
Telescope plates using the Cambridge APM. The 
plates were approximately calibrated onto i?-band 
using CCD images of star fields close to each of the 
plate centres. Initially, only the 34 objects classed 
as non-stellar in R and with 17.0 < R < 20.2 were 
considered for spectroscopy (one of these turned 
out to be a misclassificd high redshift quasar) . Six 
stellar objects were also selected later as checks on 
the APM classifier. For compact sources, or those 
with identifiable central components the probabil- 
ity of a chance misidentification is low. In the sur- 
vey region the density of objects with 17.0 < R < 
20.2 is « 0.6 arcmin~ 2 . The error on a FIRST po- 
sition is < 2 , so there is only an ss 0.2% chance 
of a misidentification. For double sources with no 
central component, (18 out of 36 objects) the prob- 
ability of a misidentification is larger. For these 
objects we followed the prescription of Lacy et al. 
(1993) by searching an ellipse with a major axis 
equal to the radio source size d and minor axis of 
d/2. The number of objects falling into the search 
region by chance, N is given in Table 1, and is 
< 1 for all our objects, and <C 1 for most of them. 
There are very unlikely to be any z < 0.7 quasars 
in the sample - the fraction of quasars in complete 
samples drops rapidly at radio luminosities below 
L 1A ~ lO^/i^WHz^sr- 1 (Willott et al. 2000). 

Spectra were obtained on the Shane 3-m Tele- 
scope at Lick Observatory on 1999 October 13-14, 

1999 November 12-14 and 1999 December 11-12 
and at the Nordic Optical Telescope (NOT) on 

2000 January 5 (all dates UT). The plate cali- 
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bration observations were made on the Shane on 
1999 September 15. The Kast spectrograph was 
used for all observations on the Shane, and the An- 
dalucia Faint Object Spectrograph (ALFOSC) for 
the NOT observations. Redshifts were determined 
from emission lines where present, otherwise ab- 
sorption features were matched both by eye and by 
a cross-correlation of the spectrum with that of a 
nearby elliptical; the dominant error sources were 
the accuracy of the wavelength calibration, and 
centroiding of spectral features, both « 0.3nm, re- 
sulting in a typical redshift error of ±0.001. Full 
details and the spectra will be presented in a fu- 
ture paper. 

3. Analysis 

To measure the two-point correlation function 
simulated surveys were constructed with the same 
selection function as the original survey, but with 
randomly assigned coordinates. The redshift se- 
lection function was estimated by using the scat- 
ter in the observed R — z relation to estimate the 
probability of an object with redshift z having a 
magnitude within the selected range (Fig. 1). A 
least deviation fit to the R — z relation for z < 0.7 
radio galaxies from the 8C-NEC redshift survey 
(Lacy et al. 1999) was used to establish a mean 
R — z relation: 

R = 21.3 -5.9 Igz . 

This is close to that derived for the 100-timcs 
radio-brighter 3C sample (Eales 1985). The scat- 
ter about this relation was measured to be 0.57 
mag, using the actual R and z values obtained for 
sources in the survey. A selection function was 
then constructed by taking the luminosity-density 
evolution model of DP and integrating to obtain a 
redshift distribution. This was then multiplied by 
the effect of the photometric selection, which was 
modelled as a pair of oppositely-tailed error func- 
tions with a width in log redshift corresponding to 
the scatter in the R — z relation, and half-power 
points corresponding to the predicted redshifts for 
objects at the bright and faint magnitude limits of 
the survey. 

In practice an upper cutoff of z = 0.45 was 
placed on the sample. Up to this redshift the 
agreement with the predicted redshift distribution 
is good, with 28.7 objects predicted in the range 
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Fig. 1. — The selection function for the sample. 
The observed redshift distribution is indicated by 
the histogram, the theoretical selection function 
by the solid line. The dashed lines indicate the 
redshift range over which we believe our complete- 
ness is good (0.19 < z < 0.45). 

0.19 < z < 0.45 compared to 29 observed. Above 
z w 0.45, however, the predicted selection func- 
tion and the observed redshift distribution diverge 
rapidly (Fig. 1), and of 47.5 objects predicted to 
be present in the redshift range 0.19 < z < 0.65, 
only 36 are found. This could be due either to the 
mistaking of high redshift galaxies for quasars by 
the APM classification program, a problem with 
the assumed luminosity function or a genuine un- 
derdense region. 

We next estimate the possible effect of incom- 
pleteness due to misclassification at z < 0.45. 
There are 12 objects within the survey selec- 
tion criteria but with stellar identifications on the 
UKST plates for which spectra have not yet been 
obtained. With one exception all have radio mor- 
phologies more consistent with being quasars than 
galaxies (i.e. unresolved, triple or core-halo struc- 
tures). Of the 13 R plate stellar identifications 
in the range 17.0 < R < 20.2 for which we have 
spectra in this region of sky, from both this work 
and the FIRST spectroscopic database, three are 
in fact galaxies, but only one is at z < 0.45. We 
therefore think it very unlikely that more than 2-3 
objects below z ~ 0.45 are missing from the sam- 
ple other than those accounted for in the selection 
function. 
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Table 1 
Radio galaxies in the sample 



Galaxy 


z 


R 


CI 


N 


010151+011414 


0.438 


19.0 


nn 


0.015 


010236-005007 


0.243 a 


17.5 


nn 


0.002 


010242-005032 


0.244 a 


17.4 


nn 


0.002 


010403-002437 


0.280 a 


17.2 


nn 


0.002 


010456+000422 


0.276 


18.4 


nn 


0.13 


010454+000400 


0.280 


17.6 


nn 


0.005 


010610+005153 


0.262 


18.2 


nn 


0.002 


010620+000843 


0.271 


18.4 


nn 


0.027 


011012-004747 


0.565 


19.6 


sn 


0.002 


0H341+010609 


0.281 


18.0 


nn 


0.002 


011425+002932 


0.355 


18.5 


nn 


0.002 


011429+000037 


0.389 


18.8 


nn 


0.017 


011527-000000 


0.381 


20.1 


nn 


0.013 


011811+010211 


0.278 


17.3 


nn 


0.002 


012012-003837 


0.236 


17.2 


nn 


0.002 


012020-002124 


0.354 


19.0 


nn 


0.002 


012030-001950 


0.352 


18.7 


nn 


0.002 


012101+005100 


0.237 


18.1 


nn 


0.033 


012156-002930 


0.437 


19.2 


nn 


0.23 


012411+004050 


0.366 


18.9 


ns 


0.002 


012927+000524 


0.392 


18.8 


nn 


0.002 


013352+011345 


0.308 


17.8 


nn 


0.002 


013506+011913 


0.358 


18.2 


nn 


0.004 


013602+001016 


0.344 


17.8 


nn 


0.007 


013807-01452 


0.640 


18.8 


ss 


0.010 


013942-000618 


0.197 a 


17.9 


ns 


0.002 


014057+001053 


0.523 


19.9 


nn 


0.002 


014213-001326 


0.526 


19.2 


nn 


0.002 


014227+001139 


0.326 


17.7 


nn 


0.015 


014251-000028 


0.272 


18.2 


nn 


0.30 


014300-000245 


0.428 


19.0 


sn 


0.43 


014317-011859 


0.520 


19.5 


nn 


0.013 


014347+004546 


0.218 


18.0 


nn 


0.002 


014714+005834 


0.638 


19.0 


nn 


0.020 


014714-012603 


0.158 


17.0 


nn 


0.55 


014752+000658 


0.448 


19.3 


nn 


0.002 



a Redshift from other FIRST follow-up 
projects 

Note. — Entries in the column labelled 'CI' 
arc 'nn' for non-stellar on both UKST plates, 
'ns' for non-stellar on the R plate and stellar on 
the B plate, 'sn' for stellar on the R plate and 
non-stcllar on the B plate, and 'ss' for stellar 
on both. N is the expected number of objects 
with 17 < R < 20.2 falling into the search 
region by chance. 
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Fig. 2. — A representation of the three dimen- 
sional distribution of galaxies in the sample. FRI 
radio sources are shown in dark grey, FRII and 
compact sources (i.e. those barely resolved or un- 
resolved by FIRST) are shown as light grey (classi- 
fications were made on the basis of radio morphol- 
ogy, and only for sources well- resolved by FIRST). 
The diameter of the spheres is 20h~ 1 comoving 
Mpc. 

In Fig. 2 we show a three dimensional represen- 
tation of the survey volume. Although this figure 
needs to be interpreted with caution as the space 
density of objects is dropping with redshift, there 
is a clump at z « 0.28 and R.A. - 0l' l 05 m , which 
seems to be comprised mostly of FRI sources. 
There is also a looser association at z w 0.35, 
and a group of three at z ~ 0.52. The sizes of 
the associations seem to be ~ 30 — 100/i _1 Mpc, 
comparable to low redshift superclusters. 

We have estimated the statistical significance of 
our detection of clustering in the complete 0.19 < 
z < 0.45 sample (median redshift w 0.3) by bin- 
ning the distribution of pair separations in 10/i _1 
Mpc bins, and comparing with the mean distri- 
bution obtained from 10000 Monte-Carlo simula- 
tions of the survey using a x 2 test- This gives the 
probability of our distribution arising by chance as 
4 x 10~ 6 (x 2 = 112 with 53 degrees of freedom). 
However, a large part of the x 2 is contributed by 
the first bin which contains two pairs with sepa- 
rations < 5/i _1 Mpc. This is close enough that the 
members of each pair could be in the same galaxy 
cluster. We therefore removed these two from the 
first bin and recalculated the x 2 statistic, this gave 
a probability of 0.006 (x 2 = 82 with 53 degrees of 
freedom) . 

To estimate the two-point correlation function, 
the numbers of data-data (DD), data-random 



(DR) and random-random (RR) pairs were mea- 
sured from the data and the simulations and the 
two-point correlation function calculated accord- 
ing to the formula espoused by Landy & Szalay 
(1993): 



£(r) = 



DD - 2DR + RR 
RR ' 



The results for the complete 0.19 < z < 0.45 
sample are shown in Fig. 3. The correlation func- 
tion was fit over the range 0-100 /i _1 Mpc, the best 
value for r was « 17/i _1 Mpc, with a range of 5- 
24/i _1 Mpc over which the probability of obtaining 
the x 2 was >5 % (assuming Poisson errors). 

As a check on the effect of possible incomplete- 
ness we added three galaxies (the largest num- 
ber we expect based on the discussion above) to 
the sample with redshifts drawn at random from 
the selection function, and with random sky posi- 
tions within the survey region. This reduced r to 
14/i _1 Mpc. Incompleteness is thus unlikely to af- 
fect our estimate by more than the random error. 
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Fig. 3. — The two-point correlation function of 
the 0.19 < z < 0.45 radio galaxies. Error bars are 
Poisson. The dotted line is the local radio galaxy 
correlation function from Peacock & Nicholson 
(1991), the dashed line the normal galaxy correla- 
tion function at z ~ 0.3 (Small et al. 1999). 



4. Conclusions and future surveys 

We have succeeded in developing an effective 
method for studying the clustering of moderate 
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redshift radio galaxies directly, and have detected 
clustering of radio galaxies at z ~ 0.3. The ampli- 
tude of the cross-correlation function we measure 
is consistent with that for radio galaxies locally. 
This is as expected in the simple model discussed 
in the introduction, in which radio source hosts 
evolve little with redshift, and is higher than that 
for normal galaxies at z ~ 0.3, for which Small et 
al. (1999) measure r = 3.7 /i _1 Mpc. At present, 
however, the small size of our survey prevents us 
ruling out all but the most extreme evolution in 
the correlation function. Expansion of this survey 
to > 100 redshifts will allow a measurement of the 
two-point correlation function to be made which 
has comparable accuracy to that for normal galax- 
ies at these redshifts. The relatively large volume 
probed by a larger survey will allow us to define a 
sample of superclusters and to examine their struc- 
tures and the evolution of those structures to the 
present day. 

The discovery of redshift clustering in deep pen- 
cil beam galaxy surveys, e.g. that of the HDF 
(Cohen et al. 2000) has raised the possibility that 
large-scale structures continue to be present in the 
Universe at least to z ~ 1. Therefore the evolu- 
tion of these structures should place interesting 
constraints on cosmology. Crucially, however, be- 
cause we can trace the evolution of radio galaxy 
hosts to z ~ 3, we can, in principle, predict how 
the bias should evolve with redshift, removing an 
important uncertainty from the interpretation of 
the results of correlation function studies. 
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